Shell filling in closed single-wall carbon nanotube quantum dots 
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We observe two-fold shell filling in the spectra of closed one- dimensional quantum dots formed 
in single-wall carbon nanotubes. Its signatures include a bimodal distribution of addition energies, 
correlations in the excitation spectra for different electron number, and alternation of the spins of 
the added electrons. This provides a contrast with quantum dots in higher dimensions, where such 
spin pairing is absent. We also see indications of an additional fourfold periodicity indicative of 
K-K' subband shells. Our results suggest that the absence of shell filling in most isolated nanotube 
dots results from disorder or nonuniformity. 

PACS numbers: 73.22.-f, 73.23.Hk, 73.63.Fg 
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Quantum dots are of broad fundamental interest, both 
because they are ubiquitous in nanoscale systems, and 
because they enable extensions and generalizations of the 
physics of atoms, molecules, nuclei, and impurities in 
solids Q. Two-dimensional (2D) dots formed in semi- 
conductor heterostructures, with controllable size and 
contact transparency, have been studied intensively for a 
decade [|J . Recently, thanks to the arrival of new systems 
including nanoparticles ||, fullerenes and nanotubes 
||, it has become possible to investigate the influences 
of geometry, band structure, atomic confirmation, vibra- 
tions, and surface chemistry on quantum dots. As with 
atoms and nuclei, in quantum dots symmetry-related or- 
bital degeneracies combined with spin and the Pauli prin- 
ciple can lead to electronic shells. Shells have indeed been 
observed in small, symmetric 3D and 2D dots, such as 
warm metal clusters || and small vertical semiconductor 
dots JtJ, respectively. However, shell filling is disrupted 
if the spatial symmetry is imperfect or if the number 
of electrons N is so large that shells overlap in energy. 
This explains why the spectra of metal nanoparticles and 
larger semiconductor dots appear chaotic, lacking (as a 
result of exchange interactions) even the two-fold shells 
that could arise from the spin degeneracy of each orbital 

One-dimensional (ID) quantum dots are now available 
|To| , [TT| , in the form of single- wall carbon nanotubes. In 
nanotubes the sole orbital symmetry is a two-fold one, 
corresponding to the K-K' subband degeneracy and re- 
sulting from the equivalence of the two atoms in the 
primitive cell of the graphene structure. Combined with 
spin this gives a possibility of four-electron shells. As yet 
unpublished reports [jl2| |l3| indicate that four-electron 
periodicity can indeed be seen in 'open' nanotube de- 
vices, where the contacts are highly transparent and 
the electron states are not localised on the tube. Here 
we report the perhaps more surprising result that even 
clean 'closed' nanotube dots, showing complete Coulomb 
blockade (CB), can exhibit simple two-electron spin shells 
in the addition spectrum, as well as the disordered reni- 
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FIG. 1: (a) Conductance vs gate voltage for a nanotube quan- 
tum dot at T = 300 mK. The conductance at room temper- 
ature was 0.3 e 2 /h. (b) Coulomb blockade peaks over a nar- 
rower range of V g . (c) Spacings AV g of the peaks in (b). The 
index n counts the added electrons relative to an arbitrary 
zero. For clarity only even n's are indicated, (d) Histogram 
of the data in (c). The sketches to the right indicate the inter- 
pretation of the two peaks in the histogram (see text). Inset: 
device structure, incorporating an atomic force microscope 
image. 



nants of four-fold shells. The results imply that, for mod- 
erate disorder, exchange corrections are very small in this 
ID system, in contrast with the 2D and 3D cases. 

Our devices were made by evaporating metal (typi- 
cally 25 nm gold on 5 nm chromium) contacts, patterned 
by electron beam lithography, on top of nanotubes grown 
by laser ablation Q and deposited from a sonicated sus- 
pension in dichloroethane onto Si02- Their geometry is 
indicated on the right of Fig. 1. The separation of the two 
contacts (source and drain) is L ~ 300 nm, and the highly 
doped silicon below acts as a metallic gate electrode. In 
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all cases the quantity measured was the dc source-drain 
current / with a bias V applied to the source and the 
drain grounded. At low temperature the majority of de- 
vices show CB oscillations of the linear conductance G as 
a function of gate voltage V g , with varying degrees of reg- 
ularity [ fl5| . We focus here on one quantum dot, selected 
from about fifty, which exhibited particularly long and 
regular sequences of CB peaks, as illustrated in Fig. la. 

A region containing 18 peaks is expanded in Fig. lb. 
The spacing, AV g , of each pair of adjacent peaks is plot- 
ted in Fig. lc. We see that it alternates, in such a way 
that we can index the blockade regions with an integer 
n which is even for all the larger AV g values. In a his- 
togram of AV 9 (Fig. Id) the odd-n values all fall in a bin 
of width 10 mV centered on 100 mV, while the even-n 
values are distributed between about 110 and 140 mV. 

The quantity AV g reflects the addition energy, ie, the 
extra energy required for adding the A^+l'th electron to 
the dot relative to the A^'th |l], [ll|. Its even-odd alter- 
nation naturally suggests two-electron shell filling. To 
analyse it further, we start with the constant-interaction 
(CI) model of a quantum dot. This model assumes that 
the total energy depends only on a given set of orbital 
energies and a single constant interaction parameter U , 
which is the repulsion between any two electrons on the 
dot. For odd N, the A^+l'th electron enters the same 
orbital as the iV'th, with energy e i; and the resulting 
separation of the CB peaks is AV g — U/ea, where e is 
the electronic charge and a < 1 is a capacitance ratio. 
For even N the A^+l'th electron enters the next orbital, 
with energy e,_|_i, and AV g = (U + ej+i — e^jea.. If we 
identify even N with even n in the experiment, then the 
interpretation of the peaks in the histogram of AV g is as 
indicated in the sketches to the right of Fig. Id. Using 
a = 0.10 for this device (which can be obtained from the 
nonlinear measurements discussed below), we then de- 
duce that U ~ 10 meV, while e^+i — e,; is distributed over 
the range 1-4 meV. These values are consistent with the 
usual results for closed nanotube dots (l7) that U/S ~ 5 
and 5 = }ivf/^L, where 5 is the mean level spacing (as- 
suming the only degeneracy is spin), vp ss 8 x 10 5 ms" 1 
is the Fermi velocity, and 5 ~ 2.7 meV for L ~ 300 nm. 

The CI model predicts that successive CB peaks should 
show related excitation spectra. The excitations can be 
studied M by making a greyscale plot of dJ/dV against 
V g and the bias V, as in Fig. 2a. In the diamond-shaped 
blockade regions no current flows and N is fixed, while 
within the cross-shaped regions the patterns reflect the 
spectra of transitions between N and N+l electrons. 
The first prediction is that for a given shell, the al- 
lowed transition energies for adding or removing electrons 
should be identical. Indeed, we observe that the crosses 
separated by odd-n diamonds, indicated by dotted lines, 
bear a clear resemblence (though they are never identi- 
cal). The second prediction is that for adjacent shells, 
the allowed transitions should follow the same pattern 
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FIG. 2: (a) A greyscale plot (darker = more positive) of 
dl/dV vs V g and V at T = 100 mK and 5=0, for the 
same device on a different cooldown. The index n counts 
electrons added relative to the leftmost diamond. The dotted 
lines have equal length, (b) The two central crosses from (a) 
are shown again, but with a relative vertical displacement. In 
the CI model the level spectrum of the dot can be read off 
from the bias at the points where the excitation lines cross 
the edge of a diamond. The resulting spectrum, given by the 
dotted horizontal lines, is almost the same for both crosses. 
The level singled out by the dashed line produces the lowest 
electron-adding transition on the right and the second lowest 
on the left, as indicated in the sketches on either side. 



but displaced such that the ground state for one corre- 
sponds to the lowest excited state for the other |18j . This 
prediction also holds true, albeit only approximately, in 
the nanotube dot. A particularly close match occurs for 
the two central crosses in Fig. 2a. They are shown again 
in Fig. 2b with the right-hand cross displaced vertically, 
allowing us to draw horizontal dotted lines which match 
all the visible electron-adding transitions on both crosses 
simultaneously. 

Furthermore, in the CI model an applied magnetic field 
B splits the two spin levels of each orbital by the Zee- 
man energy. The result is that most transitions split into 
pairs, because an electron can be added to an empty or- 
bital, or removed from a full one, with either spin. The 
only transitions which do not split are those bordering 
the odd-A^ diamonds, which involve the singly occupied 
orbital jl9| . We see this pattern whenever the peak pair- 
ing can be identified (as well as most the of the time 
even when it cannot). For example, Fig. 3 shows the 
central region from Fig. 2a at B = 6 T. Nearly all the 
visible transitions are split into pairs, with the exception 
of those bordering the n — 3 and 5 diamonds. 

The essence of the above results is captured in the 
sketches above Fig. 3, which indicate transitions corre- 
sponding to the same orbital on four consecutive CB 
peaks (crosses). Such an all-round agreement with the 
CI model has not been found in semiconductor (2D) or 
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FIG. 3: The central region of Fig. 2a in a magnetic field 
B — 6 T perpendicular to the tube axis and at T — 100 mK. 
Only the edges of the odd-n diamonds are not Zeeman split. 
This is in agreement with the constant interaction model, as 
indicated in the sketches above (see text). 

metal (3D) quantum dots. This is believed to be because 
of exchange terms, omitted from the CI model, which 
modify the level spectrum as N changes ||. In 2D and 
3D these terms are large enough (relative to S) to shuffle 
the order in which states are occupied and preclude spin 
shell filling. Our results imply that in ID dots the op- 
posite is true: exchange effects are small enough for spin 
shells to exist. 

If exchange corrections are indeed small, then the dis- 
tribution of AV g for even n in Fig. Id reflects the distri- 
bution of (^-independent) level spacings, e,+i — e,. For 
chaotic quantum dots, this distribution is given by ran- 
dom matrix theory (RMT) Q. RMT is thought to be 
appropriate for 2D and 3D dots at large N. However, it 
is not applicable for a clean ID dot, which has uniformly 
spaced levels for all N. In fact, a likely origin of the 
broad distribution we observe is suggested by scanned- 
probe measurements |pl| , which reveal that nanotubes 
contain defects. Each defect backscatters resonantly as 
a function of gate voltage in a manner dependent on its 
nature p^] . A random collection of defects can therefore 
lead to a complicated spectrum, with a corresponding 
distribution of level spacings, unique to each nanotube. 
It will however probably not randomise the Hamiltonian 
matrix sufficiently for RMT to be applicable. 

In nearly all nanotube dots studied to date (see e.g. 
refs. |l0| and Jill), most of the Coulomb diamonds are 
distorted, with accompanying complex and variable exci- 
tation patterns. Similar distortions occur in our selected 
device at some gate voltages, as can be seen in Figs. 4a 
and b. The lines of shallower slope in this grayscale plot 
correspond to a different source-gate capacitance ratio, 
implying that the geometry of the nanotube dot is not 
perfectly well defined and can effectively differ between 
eigenstates. This too can be explained by imperfections: 
a resonant defect may leave most electron states delocal- 
ized throughout the tube, but near resonance it can effec- 
tively divide the tube electronically into smaller dots with 




FIG. 4: (a) and (b) Greyscale plots of dl/dV vs V g and 
V at T = 300 mK and B = 0, for the same device on a 
different cool-down. The dotted lines indicate apparent four- 
fold grouping of the peaks, (c) Peak spacings for (a), with 
values for n = 4p plotted as solid circles, where p is an integer 
and the zero of n is chosen conveniently. 



different capacitances and very asymmetric couplings to 
the source and drain [ pl| . It may be that our device 
acts as a single quantum dot over large regions of V g 
because it contains relatively few defects. Other types 
of imperfection are also present, including a nonuniform 
(and possibly dynamically varying) electrostatic poten- 
tial along the tube, resulting from the electrode poten- 
tials and from nearby trapped charges, which it has been 
suggested Q may account for the occasional obser- 
vation of sequences of electrons appearing to enter the 
tube with the same spin p5[ . 

The observation of spin shells prompts the question 
of whether the K-K' band degeneracy also affects the 
shell structure. Combined with spin this might produce 
a fourfold, instead of a twofold, grouping of the peaks. 
No fourfold grouping is obvious in the data of Figs. 1- 
3. This does not seem surprising, because the K and K' 
states are very likely to be mixed, either by defects or at 
the contacts, lifting their degeneracy. Indeed, the mea- 
sured average level spacing only agrees with theory if this 
degeneracy is assumed lifted (leading to 5 = hvp/4L). 
However, consider instead the selections of data in Fig. 4. 
In Fig. 4a the peak spacings, plotted in Fig. 4c, are larger 
for every fourth diamond. Meanwhile, in Fig. 4b, we see 
that certain features repeat every four peaks; for exam- 
ple, the pronounced skewing of the ground-state tran- 
sition line, which we have marked with solid circles in 
Fig. 4b. The fact that a fourfold periodicity survives in 
this complex data could be explained as follows. Assume 
that a given (i'th) quartet of peaks is associated with 
wavefunctions having i nodes along the tube. The wave- 
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functions and energies may be disrupted by disorder and 
exchange. However, the i+l'th quartet will be disrupted 
in a similar way to the i'th, because the addition of one 
extra node has little consequence as long as the energy 
is not close to a defect resonance (where the scattering 
phase varies rapidly with energy). The result will be 
a complex pattern which repeats, approximately, every 
time four electrons are added. 

We have noticed a clearer four-fold peak grouping in 
devices which exhibit strong Kondo resonances of the 
conductance [ p6| . Moreover, very recent studies of similar 
devices with high contact transparency have shown un- 
mistakable four- fold patterns JlJ, |l3| . With such "open" 
contacts the current is never strongly blockaded, and the 
electronic states are only weakly confined in the tube. 
This makes them less sensitive both to disorder and to 
interactions in the tube, the latter being confirmed by 
the greatly reduced values of U/6 observed. With nearly 
ideal contacts, it appears in fact that the charging energy 
U can be neglected completely and the device treated as 
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at the contacts only p7[|. 

The electron system in a nanotube is nevertheless ex- 
pected to show strong correlations described by Luttinger 
liquid theory, for which there is a growing body of evi- 
dence [ jl7| . Many intriguing questions remain as to how 
the correlations relate to the excitation spectrum, the 
disorder, and the apparently small exchange terms that 
allow spin shells to be seen. Nanotube devices contain- 
ing fewer or no defects, and with improved control of the 
electrostatic potential, will be required to further our un- 
derstanding of these finite ID quantum systems. 

In summary, we find that closed one-dimensional quan- 
tum dots, unlike their higher dimensional counterparts, 
can exhibit twofold spin shell filling. We thank Nick 
d'Ambrumenil, Karsten Flensberg, Poul Erik Lindelof, 
Boris Muzykantskii, and Neil Wilson for helpful discus- 
sions. The experiments were supported by the Danish 
Research Councils (SNF, STVF). 
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